Resonant production of diquarks at high energy pp, ep and e + e colliders 
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Resonant productions of the first generation scalar and vector diquarks at high energy hadron- 
hadron (pp), lepton-hadron (ep) and lepton-lepton (e + e~) colliders are investigated. Taking into 
account the hadronic component of the photon, diquarks can be produced resonantly in the lepton- 
hadron and lepton-lepton collisions. Production rates, decay widths and signatures of diquarks are 
discussed using the general, SU(3)c X SU(2) w X U(1)y invariant, effective Lagrangian. The corre- 
sponding dijet backgrounds are examined in the interested invariant mass regions. The attainable 
mass limits and couplings are obtained for the diquarks that can be produced in hadron collisions 
and in resolved photon processes. It is shown that hadron collider with center of mass energy 
^J~s~ = 14 TeV will be able to discover scalar and vector diquarks with masses up to mnQ=9 TeV for 
I/"") ■ quark-diquark-quark coupling odq~0.1. Relatively, lighter diquarks can be probed at ep and e + e~ 

collisions in more clear environment. 
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Diquarks can occur in many scenarios which involve new physics beyond the standard model (SM), e.g., composite 
models Q and superstring- inspired E§ models 0. These particles can transform as anti-triplet (3*) or sextet (6) 
\ under SU(3). Diquarks carry baryon number \B\ = 2/3 and couple to a pair of quarks. They have integer spin (scalar 
or vector) and have electric charges \Q\ = 1/3, 2/3 or 4/3. 

The Collider Detector at Fermilab (CDF) has set limits on the masses of a class of scalar diquarks decaying to 
dijets with the exclusion of mass range 290 < vtluq < 420 GeV |3|, which are expected to be approximately valid for 
s ! ' other scalar diquarks. There are also indirect bounds imposed on couplings from electroweak precision data from 
LEP e + e~ collider where these bounds allow diquark-quark couplings up to a value aog ~ 0.1. 

Three types of the colliders related to the energy frontiers in particle physics research seem to be promising in 
the next decade. Namely, they are Large Hadron Collider (LHC) with the center of mass energy -y/s=14 TeV and 
luminosity L= 10 34 - 10 35 cm" 2 s _1 , International Linear Collider (ILC) with = 0.5 TeV and L= 10 34 - 10 35 
i-£h \ cm -2 s _1 , Compact Linear Collider (CLIC) with *J~s — 3 TeV and L= 10 34 — 10 35 cm~ 2 s _1 in the most preferable 
design, and the linac-ring type ep colliders, when the linear collider is constructed near the proton ring, i.e., ILCigLHC 
PL,; based ep collider with ^fs = 2.64 TeV and CLIC^LHC based ep collider with y/s = 3.74 TeV or */s = 6.48 TeV, 
having a luminosity L= 10 31 — 10 32 cm -2 s -1 . Even though the last one has a lower luminosity it can provide better 
1 conditions for investigations of a lot of phenomena comparing to ILC due to the essentially higher center of mass 
^ , energy and LHC due to more clear environment. The high energy linear e + e _ colliders have been proposed as the 
instruments that can perform precision measurements that would complement those performed at the LHC. The 
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TABLE I: The main parameters of the future ep and e + e~ colliders, L lnt denotes the integrated luminosity for one working 
year. 
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FIG. 1: Schematic presentation of resonant production of scalar or vector diquarks at three types of colliders: a) pp collisions, 
b) ep collisions with one resolved photon process and c) e + e~ collisions with two resolved photon process. 



diquarks are expected to be easily observable at LHC through their resonant production and subsequent decay to two 
jets. If relatively light diquarks are observed at LHC they are expected to be in the reach of the future linear e + e~ 
and linac-ring type ep colliders. The main parameters of future ep and e + e~ colliders are given in Table [I] 

The production and possibility of detection of diquarks have been analysed for e + e~ ||, pp [H pi and pp 0, |j| 
colliders. The single and pair production of scalar diquarks in ep collisions have been analysed in |lfjj without taking 
into account hadronic structure of the photon. Although the photon is the gauge particle of the electromagnetic 
interactions and thus pointlikc, it is known to behave like a hadron if it interacts with other hadrons. This can 
be described by the QCD-corrected quark parton model if the photon is probed at a large momentum scale. The 
importance of resolved photon contributions has been demonstrated by the DELPHI and OPAL ^1 Collaborations 
in obtaining interesting limits on leptoquark properties from ej production of leptoquarks |13| . The relevance of the 
photon substructure increases with increasing center of mass energy and becomes therefore important for forthcoming 
ep and e + e~ colliders. 

In this work, we investigate the potentials of high energy pp, ep and e + e~ colliders to search for scalar and vector 
diquarks via their resonant production subprocesses 
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where g 7 is the quark from resolved photon. The second process can also be considered as a contribution to the 
single diquark production at ep colliders. We analyze the relevant background processes in the cases of three- type of 
collisions. Schematic presentation of resonant production of diquarks at three types of colliders is shown in Fig. ^ 
(a-c). 



II. INTERACTION LAGRANGIAN 

A model independent, baryon number conserving, most general SU(3)c X SU(2)w x U{1)y invariant effective 
lagrangian for scalar and vector diquarks has the form 0, 

£|fl|=o = HlIlI^lDQ^ + (/ifldflVdfl + f' lR u R ru R )DQ% 

+ fiRUjvfd a DQ C llt + hm L T^qL • DQ| M 

+ f2q L iT 2UR.DQ2 + hq L iT2d R DQ 2 + H.c. (4) 
£|B|=2/3 = (giLql^qL + giRU c R d R )DQ1 +giRct R d R DQ 1 

_ . lc _ 

+ 9ir u rUrDQ 1 + g 3 Lq c L iT2Tq L ■ DQ3 

+ g 2 t L rdRDQ^ + hlf^URDQlv + H.c. (5) 

where the \B\ = diquarks are familiar fields, as they resemble the electroweak gauge vectors, and the neutral and 
charged Higgs scalars. Here, we consider only the \B\ = 2/3 diquarks. In Eq. (4) and (5), q^ = (ul,(Il) denotes the 
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TABLE II: Quantum numbers of the first generation, color 3 diquarks described by the effective lagrangian in the text according 
to SU(3)c x SU(2)w x U(1)y invariance where the hypercharge Y = 2(Q — 1$). 
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left handed quark spinor and q c = Cq T (q c = —q T C x ) is the charge conjugated quark field. Following H,0J, the 
possible scalar and vector diquarks are SU (3)c anti-triplets and sextets. For the sake of simplicity, color and generation 

indices are ommitted in (4) and (5). Scalar diquarks DQi, DQ l , DQ 1 are SU(2)w singlets and DQ 3 is SU(2)w 
triplet. Vector diquarks DQ 2 and DQ 2 are SU(2)w doublets. At this stage, we assume that each SM generation 
has its own diquarks and couplings in order to avoid flavour changing neutral currents (FCNC). Furthermore, the 
members of a given multiplet are assumed to be the mass degenerated. A general classification of the first generation, 
color anti-triplet (3*) diquarks is shown in Table ITU 

For the present analysis, we consider the color 3* scalar DQi or DQ® diquarks coupled to ud pairs, DQ 1 or DQ$ 

—yj + . . .1 2 

diquarks coupled to dd pair and DQ 1 or DQ^ diquarks coupled to uu pair. The vector diquarks DQ 2 and DQ 2 of 

„ 1 

type ud, DQ 2 of type dd and DQ 2 of type uu are considered. The interaction between the diquark and quark pair is 
described by the effective lagrangian (5) with different couplings. 



III. DECAY WIDTHS 



For the decay width calculation, we take the coupling as g^g — g\ + g\ for each diquark type. For numerical 
results, we will use the definition gj^g — ^ko-dq when only one type of coupling assumed to be nonzero. Diquarks 
decay into quark pairs, and the decay width T^q derived from the same lagrangian is 

r S F sg 2 D QfnDQ ,F s m DQ 

D <1 = lrfc ( ITeV } a ° Q = (6) 

r v F s9 2 DQ m DQ , T( F s m DQ 

DQ = afc 7GeV( TTeV } a ° Q = (7) 

where F$ contains the color factor for a representation including the statistical factors associated with the presence 
of identical fermions in the final state, and tudq is the mass of scalar (S) or vector (V) diquark. 



IV. SIGNAL AND BACKGROUND 



The signal for diquark production would clearly manifest itself in two jets cross sections. The differential cross 
section for any type of scalar and vector diquark resonant production can be written as follows 
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FIG. 2: Total cross sections in pp collisions for a) scalar and b) vector diquarks for different charges, with coupling strength 
cdq = 0.1, depending on their masses. 
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In the narrow width approximation {Tdq /ttidq < 0.1), the cross section of the s— channel diquark production can be 
obtained as 
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where s" is the Mandelstam variable corresponding to the square of center-of-mass energy for the subprocess. 
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A. pp Collider 

The total cross section for the resonance production of scalar diquarks at pp collider is given by 



rl dx 

r > = I —fq/p{x,Q 2 p )f q '/ P {rn 2 DQ /xs,Ql)d{s) 

'm 2 DQ /s x 



(12) 



where f q / p (x, Q 2 ) and f q i / p (x, Q 2 ) correspond to the quark distribution functions from the proton and we have used 
CTEQ5L (14j parametrization with Q 2 = s . As a consequence of this energy scan there will be quarks whose energies 
adequate for the resonance production of diquarks. The cross section is plotted against the diquark mass in Fig. [21 
(a-b) for LHC energy (y/s = 14 TeV) and coupling ang = 0.1. From these figures we find that diquarks with charge 
|Q| = 4/3 have the largest cross sections when compared to the other types. 

The scalar and vector diquarks will decay via DQ — * qiqj. Therefore, the relevant signal will be a pair of hard jets 
in the final state. At the LHC energy, major QCD background processes contributing to two-jets (2j) final states and 
their integrated cross sections are given in Table llTTI 
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TABLE III: The cross sections (in pb) for QCD backgrounds contributing to two jets final states at parton level generated by 
CompHEP with various pr cuts. 



Process 


p T > 0.1 TeV 


p T > 0.5 TeV 


p T > 1 TeV 


p T >2 TeV 


gg -> gg 


6.3xl0 5 


2.0xl0 2 


2.3x10° 


5.7xl0" 3 


qg — > qg 


6.4xl0 5 


4.8xl0 2 


l.OxlO 1 


5.7xl0~ 2 


qq -> qq 


1.0x10 s 


1.8xl0 2 


6.7x10° 


8.8xl0" 2 


gg -» qq 


2.4xl0 4 


9.8x10° 


l.OxlO -1 


2.9xl0~ 4 


qq > q'q 7 


1.6xl0 3 


2.8x10° 


1.3xl0 _1 


l.lxlO" 3 


qq — > gg 


1.5xl0 3 


2.5x10° 


6.7xl0" 2 


8.5xl0" 4 


Total 


1.4xl0 6 


8.8xl0 2 


1.9X10 1 


1.5xl0 _1 




1 23456789 
m. (TeV) 



FIG. 3: Dijet invariant mass distributions for pp — » 2jX. Resonance peaks are shown for scalar and vector diquark masses 1, 
3, 5, 7, and 9 TeV for comparison with smooth QCD backgrounds. 

The values in Table ITTT1 have been generated by CompHEP program at parton level with various pr cuts on the 
jets. It is clear that higher px cuts reduce the background cross sections significantly. These pr cuts can be translated 
into the rapidity cuts according to the relation between the pr of a jet and the rapidity with pt — rrijj/2coshy. 
Standard kinematic relations for the invariant mass rrijj and pr distributions of two-jet final states can be found in 
0- El • Figure |3| shows the dijet invariant mass distribution for the process pp — > 2j + X including the signal and the 
QCD backgrounds at the LHC. For comparison signal peaks for scalar and vector diquark masses mr>Q — 1, 3, 5, 7, 9 
TeV and oidq = 0.1 are shown on the smooth background distribution. 

In order to obtain the observability of diquarks at LHC we have calculated signal (S) and background (B) event 
estimations for an integrated luminosity of 10 5 pb -1 . The signal generated by a diquark of mass mrjQ and decay 
rate Tdq is calculated integrating the differential cross section in the two-jet invariant mass interval mryq — ^dq < 
m jj < m dq + Tdq which embraces approximately 95% of the events around the resonance. For a realistic analysis 
of the background events we take into account the finite energy resolution of the LHC- ATLAS hadronic calorimeter 
[l6l| as 5E/E — 0.5 /y/E + 0.03 for jets with \y\ < 3. The corresponding two-jet invariant mass resolution is given 
approximately by 5m jj — 0.5^/77777+ 0.02rrijj. The background is calculated by integrating the cross sections in the 
range todq — Am < rrijj < vtidq + Am with Am = max(r£>Q, 5m jj). The significance of signal over background is 
defined as S/V~B- In Fig. 01 we present S/y/~B as a function of the diquark mass for the scalar and vector diquarks 
with charges \Q\ = 4/3 and 2/3. 

If we take at least 25 signal events and S/y/B > 5 as discovery criteria, scalar (vector) diquarks with charge 
|Q| = 2/3 can be observed up to 7.5 (8.5) TeV. For the diquarks with charge \Q\ = 4/3 it is possible to cover mass 
ranges up to 9.5 (10) TeV at the LHC with L int = 10 5 pb -1 . For this luminosity, 10 7 scalar diquark events/year and 
10 8 vector diquark events/year are expected for tudq = 1 TeV. Our results show that even for much lower coupling 
constants as 10 -3 , diquarks should be seen at the LHC. 
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FIG. 4: The signal significances S/vB for diquarks as a function of diquark mass mnQ at the LHC. 

B. ep collider 

In order to calculate the total cross section for diquark production due to resolved photon process at ep colliders, 
we use the formula 



a = 
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(13) 



where / 7 / e (y) and f q >/-y{x/y, Q 2 ) corresponds to the photon energy spectrum and the probability density of 
finding a quark (qf) carrying a fraction x/y of the energy of the electrons in the initial beams, respectively. In (13) the 
integrals over the momentum fractions x and y are the mathematical representation of an energy scan performed by 
the partons coming from photons and protons. We use the photon spectrum / 7 / e (y) resulting from bremsstrahhmg 
[T^ l with the electron beam energy E e 
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We present the cross sections for scalar diquarks against their masses for coupling chdq = 0.1 and charges \Q\ = 
4/3, 1/3 and 2/3 using the bremsstrahlung process in Fig. [S]at CLIC<g>LHC based ep colliders with y/s ep — 6.48 TeV. 

In order to see the potentials of possible options of ep colliders with different center of mass energies, we plot diquark 
production cross sections depending on the ^/s ep - As seen from the Fig. El vector diquarks have larger cross sections 
than the scalar for the considered center of mass energy region. In principle, the vector and scalar type diquarks can 
be easily distinguished by the angular distribution of produced jets, the electric charge of diquarks can be determined 
by the interactions with photon. 

In the resonant production mechanism of diquarks with resolved photon process in collision ep — ► DQ+X — > 2j+X , 
most of the partcles in the final state are expected to be lost in the beam pipe, the two exceptions being the quarks 
resulting from diquark decay. As a result, the relevant signal is a px balanced coplanar pair of jets and it will contain 
two hard jets in the final state. All the relevant interactions with the diquarks are implemented into the CompHEP 
|15| program. The decay width of diquarks and the cross sections for signal and background processes are calculated 
with this program. 

In order to study contributing backgrounds we need to take into account all the diagrams contributing to two-jets 
final states. The process ep — > 2j + X has potentially large QCD background. At the ep colliders, major processes 
contributing to two-jets (2j) final states and their integrated cross sections are given in Table llVl 

The values in Table Hvl have been generated by CompHEP at parton level for transverse momentum cut of the jets 
Pt > 100 GeV. Fig. shows jet-jet invariant mass distribution for the process ep — > 2jX at CLIC<g>LHC based ep 
collider with ^/s^ = 6.48 TeV. 

In our analysis, we determine the expected significance of the signal over background for collider parameters given 
in Table ^ for ep colliders. It is convenient to collect the data in two-jet invariant mass bins since the signal is 
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FIG. 5: The total cross sections for the process ep — > DQ + X — > 2j + X at yjs cp = 6.48 TeV depending on the mass of scalar 
diquarks. 
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FIG. 6: The production cross sections a(ep — > DQ + X — > 2j + X) for scalar and vector diquarks with charge \Q\ = 4/3 and 
coupling otDQ = 0.1 depending the center of mass energy of the ep colliders. 



concentrated in a small region of the invariant mass spectrum. Therefore, we introduced the cut | rrijj — ttldq |< 25 
GeV which is efficient only for diquark masses between 0.5 and 1 TeV. For heavier diquarks rather broad resonances 
are expected and the cut | rrijj— vtidq |< 50 GeV becomes suitable for 1 TeV< itidq < 2.5 TeV. This range of the 
invariant mass will embrace approximately 95% of the signal events around the diquark resonance. For a realistic 
analysis we have estimated background events for ep colliders with the integrated luminosity of 10 3 pb _1 taking into 
account the energy resolution of the hadronic calorimeter. The statistical significances S/ \[B of diquark signal over 
background are shown in Fig. [S]for ILC&LHC and CLIC®LHC based ep colliders. From Fig. |SJ we find that scalar 
(vector) diquarks can be seen up to 1 (1.6) TeV at the ep collider with a/s = 6.48 TeV. A slightly lower limits can be 
obtained for the other ep options. 

If we take into account the integrated luminosity L = 10 3 pb _1 for the CLIC<g>LHC based ep collider with yfs — 6.48 
TeV, 44 scalar diquark events and 264 vector diquark events/year are expected for wtldq = 1 TeV and aoQ = 0.1. 
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TABLE IV: The major processes contributing to two jet (jj) final states with a cut pr > 100 GeV at the ep collisions with 
different center of mass energy y/s ep . 
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FIG. 7: The invariant mass distributions for the QCD background with the inclusion of scalar diquark signals at ep collisons. 
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FIG. 8: The signal significances S/\fB as a function of scalar and vector diquark mass uidq at different ep collider energies. 
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FIG. 9: The probability density for up (u) and down (rf) quarks within electron beam for bremsstrahlung photons. 



C. e e collider 

The cross section for the production of a pair of jets through the two resolved photon process is given by 



dx 



dyfq/e{x, Q*)fq>/e(y, Qy)v(s) 



(15) 



here f q / e (x, Q^) and f q i / e (y, Qz) correspond to the probability density of finding a quark q and q' carrying a fraction 
x and y of the energy of the electrons and positrons in the initial beams, respectively. In (15) the integrals over the 
momentum fractions x and y are the mathematical representation of an energy scan performed by the partons coming 
from electrons and positrons. The lower limits of the integrations guarantee the energy required for the production 
of quarks. As a consequence of this energy scan there will be quarks whose energies adequate for the resonance 
production of diquarks. For the probability f q / e (x,Qz) resulting from the convolution of the spectrum of photons 
f-y/e(x2) converted from the initial beams and photonic parton distribution function f q /y(xi, Q^) [l8{ . we use 
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For the relevant value entering the parton distribution in the photon, we have taken Q 7 = 



In Fig. [51 the quark probability density function f q / e (x, Q^) within electron beam are presented for bremsstrahlung 
photons at beam energy E e = 1.5 TeV. 

Here, we show only u and d quark distributions in photon since we deal with only the first family diquarks. In 
Fig. 1101 fa-b). the production cross section versus scalar and vector diquark mass is plotted for e + e~ with the center 
of mass energy ^/s e + e - = 0.5, 1 and 3 TeV. 

The two-jet background receives contributions from the annihilation process e + e" — > 7, Z — > qq and the hard 
two photon processes. The two photon processes can be of type direct, one resolved and two resolved. The last 
two processes, despite being higher order in a s , contribute significantly at lower energies. The spectrum of initial 
state radiation (ISR, photon radiation from the incoming electrons and positrons) contains the ISR scale inherent to 
the process under consideration Qjj. In our calculations, this spectrum is taken into account using the CompHEP 
program with the convolution of beamstrahlung spectra. Beamstrahlung is the process of energy loss by the incoming 
electron/positron in the field of the positron/electron bunch moving in the opposite direction. Beamstrahlung spec- 
trum is an attribute of the linear collider design and depends on the bunch geometry, bunch charge and the collision 
energy. The photon spectrum resulting from the electromagnetic interaction between the electron and positron beams 
in the intersection region is described by a more complicated expression 20] . The beamstrahlung parameters, T and 
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FIG. 10: The production cross sections for (a) scalar diquark and (b) vector diquark with charge \Q\ = 4/3 at e + e colliders 
with y/s e + e - — 0.5, 1 and 3 TeV, where the coupling qdq = 0.1. 



TABLE V: Collider parameters relevant for the calculation of beamstrahlung. T is the beamstrahlung parameter and N-, is 
the average number of photons per electron. 
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iV 7 , which in their turn are determined by the bunch design of a linear collider: 



5aNE e Ar 25a 2 N 1 

T = 5 ; r , N~ = ; (18) 

6m 3 e a z (a x + a y ) 1 12m e (a x + <x y ) Vl + T 2 /3 



where TV is the number of particles in the bunch; E e and m e are the energy and mass of the electron, respectively; 
cr x , (j y and <r z are the average sizes of the particle bunches. We presented the beamstrahlung parameters in Table 
and their effects in the background calculations in Table IVII 

We determine the expected significance of the signal over the background for collider parameters as given in Table 
[I] for ILC operating at 500 GeV, CLIC operating at 1 and 3 TeV. For the realistic analysis of the background we 
consider the generic hadron calorimeters with the two-jet mass resolution 5m jj — + 0.03mjj. In order to 

reduce the background we vetoed events exhibiting Z-bosons decaying into two-jets through the cut rrijj > 100 GeV. 
In Fig. ^2 we presented the contributions from both the annihilation process with/ without ISR+beamstrahlung and 
two-photon processes to dijet cross sections at ^/s e + e - = 3 TeV. 



TABLE VI: Background calculation with the ISR and beamstrahlung effects. The numbers are the cross sections in pb for the 
process e+e~ -> 7, Z -> qq with p 3 T > 20 GeV (p£, > 100 GeV). 



Collider 


ILC 


CLIC 




parameter 


500 GeV 


1 TeV 


3 TeV 


no ISR 


2.22(1.96) 


0.546(0.672) 


0.0604(0.0602) 


with ISR 


8.16(2.11) 


1.96(0.629) 


0.276(0.0829) 


with ISR+beamstrahlung 


8.38(2.16) 


4.30(1.04) 


67.55(1.39) 
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FIG. 11: The transverse momentum distributions of the jet. 
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FIG. 12: The scalar and vector diqurk signal significance as a function of their masses. 



We conclude that signal significances scale with the collider luminosity as yL, and higher luminosities are desirable 
for a good signal over background discrimination. In determining the background events we considered some cuts on 
transverse momentum pr > 100 GeV of the jets. In Fig. ^1 we present the signal significance that can be reached at 
e+e - collider with at ^/s e + e - = 3 TeV. 

We take the indication of signal S/ \f~B > 3 and find that vector diquarks can be probed up to 1 TeV at e + e~ 
collider with y?s e + e - — 3 TeV. For the e + e~ colliders running at the center of mass energy of 1 TeV we can search 
for scalar diquarks up to 325 GeV and vector diquarks up to 375 GeV. 



V. CONCLUSION 



The resonant production of scalar and vector diquarks at LHC have large cross section. With reasonable cuts, it 
may be possible to cover mass ranges up to 10 TeV for coupling ajj =0.1. For smaller couplings as ajj = 10 -3 , it is 
still possible to probe diquarks up to the mass of 4 TeV at an integrated luminosity L — 10 2 fb _1 . We find that vector 
diquarks can be seen up to 1.6 TeV at the ep collider with y/s = 6.48 TeV and a slightly lower limits can be obtained 
for scalar diquarks. The two hard photon processes prevail at energies below \fs/2 while the annihilation prevails at 
higher energies. Therefore, the attainable mass limits for diquarks are smaller than \/s/2 at linear colliders with the 
parameters given in Table [I] 

A limited number of measurements of diquark properties can be carried out at the hadron colliders. The future 
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high energy linear e + e~ colliders and linac-ring type ep colliders based on the hadron and lepton colliders would 
complement the measurements performed at hadron colliders. The resolved photon contributions to the search for 
new physics becomes important for forthcoming ep and e + e~ collisions. 
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